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Abstract 

Aims. The main goal of this work is to study element ratios that are important for the formation of planets of different masses. 
Methods. We study potential correlations between the existence of planetary companions and the relative elemental abundances of 
their host stars. We use a large sample of FGK-type dwarf stars for which precise Mg, Si, and Fe abundances have been derived using 
HARPS high-resolution and high-quality data. 

Results. A first analysis of the data suggests that low-mass planet host stars show higher [Mg/Si] ratios, while giant planet hosts 
present [Mg/Si] that is lower than field stars. However, we found that the [Mg/Si] ratio significantly depends on metallicity through 
Galactic chemical evolution. After removing the Galactic evolution trend only the difference in the [Mg/Si] elemental ratio between 
low-mass planet hosts and non-hosts was present in a significant way. These results suggests that low-mass planets are more prevalent 
around stars with high [Mg/Si]. 

Conclusions. Our results demonstrate the importance of Galactic chemical evolution and indicate that it may play an important role 
in the planetary internal structure and composition. The results also show that abundance ratios may be a very relevant issue for our 
understanding of planet formation and evolution. 
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1. Introduction 

It is well known that stellar metallicity is a key parameter for 
the formation and evolution of planets. In particular, it has been 
shown that giant planet formation is more efficient around metal- 
rich stars (e.g. Gonzalez|1997t Santos et al.|2004[ Johnson et al. 


|2010[ Mortier et al.||2013 l, while this planet-metallicity corre¬ 

lation probably does not hold for lower-mass/small-size plan¬ 
ets (e.g. Sousa et al.|201 lb] Mayor et al.|201 1 Buchhave et al. 
|2012| l. It has also been recently shown that stellar metallicity 
plays an important role in the architecture of planetary orbits 


(e.g. Dawson & Murray-ClaypoTSj Beauge & Nesvomy|2013[ 

[Adibekyan et al.|2013b[). ' ' 

In all of the works cited, the iron content was used as a 
proxy for the overall metallicity. Building on this, many au¬ 
thors searched for chemical peculiarities of planet hosting stars 
in terms of abundances of individual elements. While many 
contradictory results were obtained (e.g. [Santos et al.||2000^ 


Bodaghee et al.||2003[ [Bond et al.||2006[ Robinson et al. 2006 


Kang et al. 201 1[ Brugamyer et al.|2011| l, the enhancement 
of (T-elements of iron-poor planet hosts was shown to be ro- 
bust ([Haywood||2008 1 Adibekyan et al.||2012a]b i. Interestingly, 
Adibekyan et al. ^ 2012a[ l showed that low-mass/small-size plan¬ 
ets show a-enhancement at iron-poor regimes. This result sug¬ 


gests that although low-mass/small-size planets can be formed at 
a wide range of metallicities ([Buchhave et al.[[2012| [Mordasini 
et al.|2012| i, extra heavy elements are required to “compensate” 
for the lack of iron in low-metallicity environments. 


Interestingly, observational and theoretical studies suggest 
that abundances ratios of individual elements in a circumstel- 
lar disk determine the structure and composition of the plan¬ 
ets that are formed (e.g. Grasset et al. [[20091 [Bond et al. 2010 


Delgado Mena et al. 2010[ Rogers & Seager 2010 Thiabaud 
et al.|2014[[2015[[Dorn et al.|2015| l. Among the most important 
mineralogical ratios are the Mg/Si, Fe/Si, and C/O ratios that 
can be used to constrain the chemical composition of terrestrial 
planets (e.g. [Bond et al.|2010[ [Thiabaud et al.|2015[['Dorn et aL| 
20151. Recently, Dressing et al. (2015 1 showed that all the ter¬ 
restrial planets in their sample with masses lower than bM® lie 
in the same mass-radius area as Earth and Venus. Their sugges¬ 
tion that these planets can be explained by Earth-like composi¬ 
tion was successfully tested by Santos et al. ( [2015 1. The authors 
derived iron mass fractions of the planets using chemical abun¬ 
dance proxies of their host stars and found that it varies from 
27.5 to 34.7 %, being similar with the iron mass fraction for the 
Earth (29-32% e.g. [McDonough & Sun|1995l ). 


1 


















































































V. Adibekyan et al.: From stellar to planetary composition: Galactic chemical evolution of Mg/Si mineralogical ratio 




G [Mg/Si] 


o 




[Fe/H] 


Figure 1. Left. [Mg/Si] versus [Fe/H] for stars without detected planets (black dots), stars hosting low-mass planets (blue squares), 
and stars hosting high-mass planets (red filled circles). Cumulative distribution of [Mg/Si] {top right), [Mg/Si] after correcting for 
the Galactic chemical evolution {middle right), and [Fe/H] {bottom right) for stars with and without detected planets. The black 
solid line provides linear fit to all the data points and the green dashed line shows separation between Galactic thin- and thick-disk 
stars. 


Motivated by these recent works, we decided to explore these 
important elemental ratios in stars harbouring different type of 
planets. The main goal of this work is to verify whether planet 
hosts exhibit peculiar elemental ratios of Mg/Si and to under¬ 
stand the origin of possible peculiarities. We organized this 
Letter as follows. In Sect.j^ we present the sample; in Sect.[^ 
we present and discuss the results. We finish the paper with con¬ 
cluding remarks in Sect.j^ 


2. The sample 

Our initial sample consists of 1111 FGK-type stars with and 
without detected planets that were observed with the HARPS 
spectrograph ( [Adibekyan et al.|2012c| ). This work provides pre¬ 
cise chemical abundances of 12 refractory elements derived in a 
homogeneous manner. The stellar atmospheric parameters of the 
stars were taken from jSousa et al.| ( |2008|[201 lb|a| l. 

To study stars with the most precise stellar parameters and 
chemical abundances we selected only dwarf (logg > 4) stars 
with effective temperatures within 500 K from that of the Sun 


(Tq = 5777 K) (see the discussion in Adibekyan et al.||2012ct 
Tsantaki et al.|2013 for details). We also constrained our sample 


by selecting stars for which errors on Mg, Si, and Fe abundances 
are smaller than 0.2 dex. We note that this value is very con¬ 
servative and guarantees exclusion of any possible outliers for 
which very precise abundances could not be derived (due to e.g. 
cosmic rays around some spectral lines or low signal-to-noise ra¬ 
tio in the spectra). We also note that our results are not sensitive 
to this constraint on precision and all the conclusions remain the 
same. Finally, we only selected stars that have metallicity higher 
than -0.6 dex which is the lower [Fe/H] limit of planet hosting 
stars in the sample. 


Our final sample consists of 493 stars without detected plan¬ 
etary companions, 19 stars hosting super-Earths or Neptune-like 
planets (the most massive planet in the system has M^ < 3QMq) 
and 77 Jovian hosts (M^ > 3 OM 0 ). 
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Table 1. KS and AD probabilities that planet hosts and stars 
without detected planets come from the same parent distribution. 




[Mg/Si] 

[Mg/Si]„„ 

[Fe/H] 


Neptune - Non-host 

0.037 

0.009 

0.973 

^KS 

Jovian - Non-host 

0.051 

0.736 

1.2x10^^ 


Neptune - Non-host 

0.020 

0.003 

0.883 

^AD 

Jovian - Non-host 

0.017 

0.841 

8.6x10^^ 


3. Results and discussion 


Since [Adibekyan et al.| ( |2012a|c| l already discussed the possible 
differences in individual elemental abundances and [X/Fe] abun¬ 
dances ratios between stars without detected planets and stars 
hosting low-mass and high-mass planets, in this work we fo¬ 
cus on the [Mg/Si] ratio. We remind the reader that [Adibekyan 
et al.| ( |2012a l showed that stars hosting planets show enhance¬ 
ment in a-elements (high [Si/Fe] and [Mg/Fe] ratios) relative to 
non-hosts at low-iron regime. 

In the left panel of Fig. [T we plot the dependence of the 
[Mg/Si] ratio on stellar metallicity for stars with and without 
detected planets. The top right panel shows the cumulative dis¬ 
tributions of [Mg/Si] ratios for the three groups of stars. As 
can be easily seen, the [Mg/Si] distributions differ among the 
groups: Low-mass planet hosts have higher [Mg/Si] and high- 
mass planet hosts have lower [Mg/Si] ratios than their non-host 
counterparts. To quantify these differences we applied the two- 
sample Kolmogorov-Smirnov (KS) and the Anderson-Darling 
(AD|^ tests. Both KS and AD statistics yield a low probability 
that low-mass and high-mass planet hosts come from the same 
[Mg/Si] distribution as the non-host stars (see Table[2l. 

The left plot, however, clearly shows that the [Mg/Si] ra¬ 
tio depends on metallicitjj^ Taking also into account the fact 
that stars hosting high-mass planets show a metallicity distribu¬ 
tion different from that of “single” stars (see the last column of 
Table[^, one can argue that this metallicity difference can be 
responsible for the difference in [Mg/Si] distribution between 
Jovian hosts and non-hosts. To remove the trend of [Mg/Si] with 
Galactic chemical evolution (GCE) we fitted all our data points 
(all stars with and without detected planets) with a linear depen¬ 
dence on metallicity and then subtracted the fit. To the subtracted 
data we applied the KS and AD tests again, and found that the 
difference in [Mg/Si] between Jovian hosts and the single stars 
disappears, suggesting that the [Mg/Si] difference in the original 
dataset was just a reflection of GCE and the shift of giant-planet 
hosts towards high metallicities. 

The difference in [Mg/Si] after correction for the GCE 
([Mg/Sileorr) between Neptune hosts and non-hosts remains, and 
the applied statistical tests predict even lower probabilities that 
the two samples come from the same parent distribution. To 
evaluate the statistical significance of the observed difference in 
[Mg/Si]corr, we performed two simple Monte Carlo (MC) tests. 
In the first test, we randomly draw 19 stars (this is the total num¬ 
ber of low-mass planet hosts) from the non-host sample and ap¬ 
plied the KS test on the distributions. We repeated the entire pro¬ 
cess 10^ times and counted the number of trials in which the P/fj 
is equal to or smaller than 0.009 (the probability that Neptunian 
hosts and non-hosts come from the same [Mg/Si]eorr distribu- 

* AD is similar to the KS test, but more sensitive (gives more weight) 
to the tails of distribution. 

^ The performed liner fit suggests that [Mg/Si] = 0.0045(±0.0017) - 
0.0873(±0.0071)x[Fe/H]. 


tion; see Table[T]l. Only in 0.5% of the cases did we obtain a 
Pats < 0.009, which means that the observed difference is most 
likely not obtained by chance. 

We performed a second MC test in which we counted the 
number of low-mass planet hosts that lie above the linear fit of 
the GCE (15 out of 19), and again 10^ times randomly drew 19 
stars from the non-host sample every time counting the num¬ 
ber of trials in which at least 15 stars lie above this line. The 
test again showed that the probability of obtaining the observed 
[Mg/Si]corr distribution by chance is very low, i.e. 0.9%. 

A word of caution should be added at this point. The 
abundances of both Mg and Si are derived by assuming lo¬ 
cal thermodynamic equilibrium (LTE). The non-LTE effects de¬ 
pend on the stellar parameters of the stars. Hence, if the sam¬ 
ples of the Neptunian hosts and stars without detected plan¬ 
ets have significantly different distribution of stellar parame¬ 
ters, then the obtained difference in the [Mg/Si] ratio might be 
affected. Several works in the literature have already tried to 
quantify the non-LTE effects in Mg and Si, and their depen¬ 
dence on the stellar parameters (e.g. Zhao & Gehren|2000[ [Shi 


et al.[[2011t Merl e et al.|20TT| [Sukhorukov & Shchukina||26l2t 


Bergemann et al. 2013| l. These works show that the non-LTE 
effects are usually stronger, being non-negligible for evolved, 
hot, and metal-poor stars. As mentioned in Sect.j^ the stars 
in this work were selected specifically to be main sequence 
(dwarfs), and to have temperatures similar to that of the Sun. 
When comparing the atmospheric parameters of the two samples 
we found that the hosts of low-mass planets are slightly cooler 
ihTeff = (5645+196)^,^,- (5771+247),= -126 K) 
than their counterparts without planets, but have similar metal- 
licities (A[Ee/H] = (-Qm±0.l9)Neptune - (-0.08±0.23)„o„_ao,, = 
0.0 dex) and are at the same evolutionary stage (Aloggr = 
{A32>±Q.01)Neptune - {AA\±0.\2)non-host = -0.03 dex). The re¬ 


sults of the aforementioned references leads us to conclude that 
this small difference of 100 K (which is within the dispersion 
of the Tg// values of each sub-sample) in temperature is not ex¬ 
pected to be responsible for the significant difference observed in 
[Mg/Si] between Neptunian hosts and non-hosts. However, only 
the non-LTE corrections of the Mg and Si abundances would 
provide precise quantification of the impact of this difference in 
Tg/j on the [Mg/Si] ratio. 


4. Concluding remarks 

The dependence of [Mg/Si] on metallicity, i.e. the Galactic evo¬ 
lution of [Mg/Si] ratio, may play an important role in the internal 
structure and composition of terrestrial planets. Because they are 
a-capture elements. Mg and Si are mostly products of massive 
stars. The main mechanisms for the production of Mg are C and 
Ne burning, while for Si, it is explosive and non-explosive O 
burning (e.g. Woosley & Weaver] 1995 [ Thielemann et al.|1996| l. 
While Mg (like O) is almost exclusively produced in massive 
stars, the production of Si by Supernovae type la is not negli¬ 
gible. This difference in the production sites is reflected in the 
trends observed in the [Mg/Fe] and [Si/Fe] versus [Fe/H] plots 
(e.g. [Fulbright et al.|2007 Adibekyan et al.|2012c| l. 

We find that low-mass planets are more prevalent around 
stars with high [Mg/Si] and that the [Mg/Si] ratio depends on 
stellar metallicity may have important implications. The first 
finding is that Mg/Si ratio probably plays a very important role 
in the formation of these low-mass planets. The second result 
is that there is probably a dependence of the planetary struc¬ 
ture on the Galactic evolution i.e. low-mass planets that were 
formed at different times and at different places in the Galaxy 
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may tend to have different structures and compositions. In par¬ 
ticular, the left panel of Fig. [T] shows that stars belonging to the 
Galactic thick disk (the stars above the green dashed line) have 
higher [Mg/Si] ratio than their thin disk counterparts (below 
the the line) Taking into account the fact that planet-hosting 
stars at low metallicities mostly belong to the Galactic thick 
disk ( |Haywood||2008[ [Adibekyan et al.|[2012b| l one can specu¬ 
late that the oldest terrestrial planetary systems (e.g. Kepler-444; 
Campante et al.|[2015|l may have different compositions to the 


terrestrial solar system planets; however, for the exact derivation 
of the chemical properties of these planets precise abundances 
of C and O are also needecQ While we already derived oxy¬ 


gen abundance for most of the stars in this sample (Bertran de 
[Lis et al.l|2015| ) the derivation of carbon abundances is still in 
progress (Suarez-Andres et al., in prep; Delgado Mean et al., in 
prep). 

Summarizing our results, we can conclude that stars formed 
at different times and places in the Galaxy have a different prob¬ 
ability of forming low-mass planets, and the composition of the 
formed planets will also depend on the chemical composition of 
the environment in which they formed. 
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